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ABSTRACT 

We  observed  v  =  \,J  =  1-0  43  GHz  SiO  maser  emission  toward  the  Mira  variable  IK  Tauri  (IK  Tau)  using  the 
Very  Long  Baseline  Array  (VLB A).  The  images  resulting  from  these  observations  show  that  SiO  masers  form  a 
highly  elliptical  ring  of  emission  approximately  58x32  mas  with  an  axial  ratio  of  1.8:1.  The  major  axis  of  this 
elliptical  distribution  is  oriented  at  a  position  angle  of  ^59°.  The  line-of-sight  velocity  structure  of  the  SiO  masers 
has  an  apparent  axis  of  symmetry  consistent  with  the  elongation  axis  of  the  maser  distribution.  Relative  to  the 
assumed  stellar  velocity  of  35  km  s~,  the  blue-  and  redshifted  masers  were  found  to  lie  to  the  northwest  and 
southeast  of  this  symmetry  axis,  respectively.  This  velocity  structure  suggests  a  northwest-southeast  (NW-SE) 
rotation  of  the  SiO  maser  shell  with  an  equatorial  velocity,  which  we  determine  to  be  ~3.6  km  s  1 .  Such  a  NW-SE 
rotation  is  in  agreement  with  a  circumstellar  envelope  geometry  invoked  to  explain  previous  H20  and  OH  maser 
observations.  In  this  geometry,  FLO  and  OH  masers  are  preferentially  created  in  a  region  of  enhanced  density  along 
the  NE-SW  equator  orthogonal  to  the  rotation/polar  axis  suggested  by  the  SiO  maser  velocities. 

Subject  headings:  circumstellar  matter  —  masers  —  radio  lines:  stars  —  stars:  AGB  and  post-AGB  — 
stars:  individual  (IK  Tauri) 


1.  INTRODUCTION 

The  star  IK  Tau,  also  know  as  NML  Tau  ( Neugebauer  et  al. 
1965),  is  an  extremely  red  Mira  variable  with  an  unusually  late 
spectral  type  ranging  from  M8. 1  to  Ml  1 .2  (Wing  &  Lockwood 
1973).  The  period  of  IK  Tau  is  ~470  days,  as  determined  from 
infrared  photometry  measurements  at  1  jon  (Wing  &  Lockwood 
1973)  and  at  11  gm  (Hale  et  al.  1997).  Olofsson  et  al.  (1998) 
computed  a  distance  of  250  pc  from  integrated  visual,  near- 
infrared,  and  IRAS  data  using  the  period-luminosity  relation¬ 
ship  of  Whitelock  et  al.  (1994).  This  distance  is  consistent  with 
the  270  pc  determined  by  Hyland  et  al.  (1972)  using  infrared 
flux  measurements  and  an  assumed  long-period  variable  max¬ 
imum  luminosity  of  10~4.  Mass-loss  rates  estimated  from  the 
CO(J  =  1-0)  line  range  from  5.1  xl0~6  MQ  yr_1  (Knapp  & 
Morris  1985)  to  2.4x  10~6  MQ  yr~'  (Olofsson  et  al.  1998). 

Long-period  variables  such  as  IK  Tau  often  exhibit  circum¬ 
stellar  maser  emission  in  one  or  more  molecular  species  (i.e., 
OH,  H20,  or  SiO).  IK  Tau  is  one  such  star  that  has  maser 
emission  in  all  three  species.  The  OH  and  H20  masers  toward 
IK  Tau  have  been  mapped  using  connected  element  and  radio- 
linked  interferometers  (Lane  et  al.  1987;  Bowers  et  al.  1989, 
1993;  Yates  &  Cohen  1994;  Bains  et  al.  2003),  revealing  the 
structure  of  the  circumstellar  envelope  (CSE)  on  both  large 
scales  (2" -4"  traced  by  the  OH  masers)  and  intermediate  scales 
(200-350  mas  traced  by  the  H20  masers).  Such  studies  have 
been  useful  in  determining  the  structure  of  the  CSE  on  these 
scales.  The  SiO  masers,  which  lie  close  to  the  star,  provide  a 
unique  probe  of  the  morphology  of  the  CSE  and  the  velocity  of 
the  gas  in  a  region  just  a  few  stellar  radii  from  the  photosphere. 
Previous  very  long  baseline  interferometry  (VLBI)  observations 
of  long-period  variables  have  demonstrated  the  various  morphol¬ 
ogies  that  exist  in  the  SiO  maser  region.  Some  stars  exhibit  ring¬ 
like  SiO  maser  distributions  (Diamond  et  al.  1994;  Diamond  & 
Kemball  1998;  Cotton  et  al.  2004),  while  others  (e.g.,  NML  Cyg) 


show  more  elliptical  distributions  ( Boboltz  &  Marvel  2000).  Pre¬ 
vious  SiO  maser  studies  have  also  revealed  the  velocity  structure 
of  the  nearby  CSE  to  be  complex,  showing  elements  of  con¬ 
traction  (Boboltz  et  al.  1997),  expansion  (  Diamond  &  Kemball 
2003),  and  even  rotation  (Boboltz  &  Marvel  2000;  Hollis  et  al. 
2001;  Sanchez  Contreras  et  al.  2002;  Cotton  et  al.  2004).  In  this 
article,  we  present  the  results  from  our  VLBI  observations  of  the 
SiO  maser  emission  toward  the  long-period  variable  IK  Tau.  We 
discuss  how  these  observations  fit  into  the  overall  picture  of  the 
CSE  as  traced  by  the  various  maser  species  and  the  possibility  of 
axial  symmetry  and  rotation  in  the  envelope  of  IK  Tau. 

2.  OBSERVATIONS  AND  REDUCTION 

We  observed  the  v  =  \,J  =  1-0  SiO  maser  emission  toward 
IK  Tau  (a  =  03h53m28!8,  6  =  11°24'22"4  [J2000J)  on  1996 
April  5  from  20:00  to  24:00  UT.  Based  on  a  period  of 470  days 
(Wing  &  Lockwood  1973;  Hale  et  al.  1997)  and  a  maximum 
at  JD  2,439,440  (Wing  &  Lockwood  1973),  the  stellar  phase  of 
IK  Tau  was  ~0.85  during  the  epoch  of  our  observations.  IK  Tau 
and  a  continuum  calibrator  (0423-013)  were  observed  using  the 
1 0  stations  of  the  VLB  A  plus  a  single  antenna  of  the  Very  Large 
Array  (VLA).  Both  the  VLB  A  and  VLA  are  operated  by  the 
National  Radio  Astronomy  Observatory  (NRAO). 1  A  reference 
frequency  of  43. 122027  GHz  was  used  for  the  v  =  1  ,J=  1-0 
SiO  transition.  Data  were  recorded  in  dual  circular  polarization 
using  a  single  4  MHz  (27.8  km  s~*)  band  centered  on  the  local 
standard  of  rest  (LSR)  velocity  of  34.0  km  s_1.  System  tem¬ 
peratures  and  point  source  sensitivities  were  on  the  order  of 
~150  K  and  ~11  Jy  K_1,  respectively. 


1  The  National  Radio  Astronomy  Observatory  is  a  facility  of  the  National 
Science  Foundation  operated  under  cooperative  agreement  by  Associated  Uni¬ 
versities,  Inc. 
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Fig.  1. — Total  intensity  VLBI  images  of  the  v—  1,  J  —  1-0  SiO  maser 
emission  toward  IK  Tau  integrated  over  the  LSR  velocity  range  from  +28.8  to 
+41.6  km  s"1.  Contour  levels  are  -3%,  -1.5%,  1.5%,  3%,  5%,  10%,  20%, 
40%,  and  80%  of  the  peak  integrated  flux  density  of  49.9  Jy  beam-1.  The 
synthesized  beam  is  0.54x0.41  mas  at  a  position  angle  of  14?1.  The  dashed 
circle  represents  a  stellar  uniform  disk  diameter  of  20.2  mas  from  Monnier  et  al. 
(2004).  The  location  of  the  stellar  disk  is  conjecture. 

The  data  were  correlated  at  the  VLBA  correlator  operated  by 
the  NRAO  in  Socorro,  New  Mexico.  Auto-  and  cross-correlation 
spectra  consisting  of  128  channels  with  channel  spacings  of 
31.25  kHz  (~0.2  km  s  ')  were  produced  by  the  correlator. 
Calibration  was  performed  using  the  Astronomical  Image  Pro¬ 
cessing  System  (A1PS)  maintained  by  NRAO.  Total  intensity 
(Stokes  I)  data  were  calibrated  in  accordance  with  the  proce¬ 
dures  outlined  in  Diamond  (1989).  The  bandpass  response  was 
determined  from  scans  on  the  continuum  calibrator  and  used  to 
correct  the  target  source  data.  The  time-dependent  gains  of  all 
antennas  relative  to  a  reference  antenna  were  determined  by  fit¬ 
ting  a  total-power  spectrum  (from  the  reference  antenna  with 
the  target  source  at  a  high  elevation)  to  the  total-power  spectrum 
of  each  antenna.  The  absolute  flux  density  scale  was  established 
by  scaling  these  gains  by  the  system  temperature  and  gain  of  the 
reference  antenna.  Errors  in  the  gain  and  pointing  of  the  refer¬ 
ence  antenna  and  the  atmospheric  opacity  contribute  to  the  error 
in  the  absolute  amplitude  calibration,  which  is  accurate  to  about 
1 5%-20%.  Residual  group  delays  were  determined  by  perform¬ 
ing  a  fringe  fit  to  the  continuum  calibrator  sources  and  were 
applied  to  the  spectral  line  data.  Residual  fringe-rate  solutions 
were  obtained  by  fringe  fitting  a  strong  reference  feature  in  the 
channel  at  FLsr  =  36.2  km  s  1  and  were  applied  to  all  channels 
in  the  spectrum.  An  iterative  self-calibration  and  imaging  pro¬ 
cedure  was  performed  to  map  the  reference  channel  at  36.2  km 
s- 1 .  The  resulting  phase  solutions  were  applied  to  all  channels 
in  the  band. 

Images  of  the  SiO  maser  emission  consisting  of  1024  x 
1024  pixels  (~72x72  mas)  were  generated  using  a  synthesized 
beam  of  0.54x0.41  mas  in  natural  weighting.  Images  were  pro¬ 
duced  for  spectral  channels  from  28.8  to  41.6  km  s-1,  forming 
an  image  cube  of  60  planes.  Off-source  rms  noise  estimates  in 
the  spectral  channel  images  range  from  10  mJy  beam- 1  in  planes 
with  weak  or  no  maser  emission  to  80  mJy  beam-1  in  dynamic- 


Fig.  2. — Maser  spot  map  as  derived  from  the  two-dimensional  Gaussian  fits 
to  masers  in  each  plane  in  the  image  cube.  Point  sizes  are  proportional  to  the 
logarithm  of  the  flux  density.  The  ellipse  indicates  a  flux-density- weighted  least- 
squares  fit  to  the  maser  distribution. 

range-limited  planes  containing  strong  (^15  Jy  beam-1)  ma¬ 
ser  emission.  Figure  1  shows  the  total  intensity  contour  map  of 
the  43 . 1  GHz  SiO  masers  integrated  over  the  LSR  velocity  range 
above.  The  dashed  circle  in  the  figure  represents  the  stellar  photo¬ 
sphere  with  a  diameter  of  20.2  mas.  This  diameter  was  determined 
from  a  uniform  disk  ( UD)  diameter  fit  to  combined  Keck  aperture- 
masking  data  and  Infrared  Optical  Telescope  Array  ( IOTA)  data 
in  the  near- infrared  (Monnier  et  al.  2004).  It  is  apparent  from 
Figure  1  that  the  SiO  masers  toward  IK  Tau  form  an  elongated  ring 
roughly  58x32  mas  across  oriented  in  the  northeast-southwest 
(NE-SW)  direction.  If  the  star  is  assumed  to  be  in  the  center  of 
the  distribution,  maser  distances  range  from  ~1  .6  to  2.8  stellar 
radii. 

In  order  to  identify  and  parameterize  maser  components,  two- 
dimensional  Gaussian  functions  were  fitted  to  the  emission  in 
each  spectral  (velocity)  channel  using  the  A1PS  task  SAD.  Image 
quality  was  assessed  using  the  off-source  rms  noise  and  the  deep¬ 
est  negative  pixel  in  the  image.  A  cutoff,  above  which  emission 
features  were  fit,  was  conservatively  set  to  the  greater  of  10  <rrms 
or  80%  of  the  absolute  value  of  the  deepest  negative  pixel  in  the 
plane.  Errors  reported  by  SAD  were  based  on  the  source  size 
divided  by  twice  the  signal-to-noise  ratio  (S/N )  in  the  image  and 
ranged  from  5  fias  for  features  with  high  S/N  to  37 1  /xas  for  fea¬ 
tures  with  lower  S/N.  The  emission  features  identified  by  SAD 
are  represented  by  the  circles  in  Figure  2.  In  the  figure,  point 
sizes  are  proportional  to  the  logarithm  of  the  fitted  flux  density, 
which  ranged  from  0.2  to  18.3  Jy.  Comparing  Figures  1  and  2, 
we  see  that  the  identified  features  accurately  represent  the  emis¬ 
sion  summed  over  all  velocity  channels  in  the  image  cube. 

3.  RESULTS  AND  DISCUSSION 
3.1.  Axial  Symmetry  in  the  CSE  of  IK  Tan 

Figures  1  and  2  show  that  the  distribution  of  the  43.1  GHz 
SiO  maser  emission  toward  IK  Tau  was  clearly  not  circularly 
symmetric  during  the  epoch  of  our  observations.  Instead,  the 
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distribution  appears  to  be  elongated  along  an  axis  running  from 
the  northeast  to  the  southwest  in  the  images.  We  characterized 
this  morphology  by  performing  a  least-squares  fit  of  an  ellipse 
to  the  distribution  of  masers  weighted  by  the  flux  density  of 
each  component.  The  resulting  ellipse  is  plotted  in  Figure  2.  The 
lengths  of  the  semimajor  and  semiminor  axes  were  found  to  be  29 
and  16  mas,  respectively,  with  the  semimajor  axis  of  the  ellipse 
oriented  at  59°  measured  east  of  north.  At  the  assumed  distance 
to  IK  Tau  of  250  pc,  the  distribution  of  SiO  masers  is  approxi¬ 
mately  15x8  AU.  The  SiO  maser  distributions  of  5120  late-type 
stars  have  previously  been  imaged  with  angular  resolutions 
provided  by  VLBI  (e.g.,  Cotton  et  al.  2004;  Diamond  &  Kemball 
2003;  Sanchez  Contreras  et  al.  2002;  Hollis  et  al.  2001;  Boboltz 
&  Marvel  2000;  Diamond  et  al.  1994).  Of  these  stars,  IK  Tau, 
with  an  axial  ratio  of  ~1.8:1,  appears  to  have  one  of  the  most 
elongated  SiO  rings  yet  observed. 

As  mentioned  previously,  Monnier  et  al.  (2004)  were  able 
to  separate  the  near-infrared  stellar  and  dust  contributions  and 
determine  a  UD  diameter  of  20.2  mas  for  the  photosphere  of 
IK  Tau.  However,  they  found  that  the  longest  (27  m)  IOTA  base¬ 
line  data  were  inconsistent  with  this  fit,  and  the  data  were  thus 
ignored.  Possible  explanations  for  this  deviation  are  departures 
from  uniform  photospheric  brightness  (e.g.,  stellar  hot  spots  or 
extended  molecular  emission)  or  the  presence  of  a  binary  com¬ 
panion  (Monnier  et  al.  2004). 

The  circumstellar  dust  shell,  typically  outside  the  SiO  maser 
region,  has  been  observed  for  IK  Tau  using  long-baseline  inter¬ 
ferometers  operating  at  mid-infrared  (Danchi  et  al.  1994;  Hale 
et  al.  1997)  and  near-infrared  ( Monnier  et  al.  2004)  wavelengths. 
Hale  et  al.  ( 1 997 )  found  that  their  Infrared  Spatial  Interferometer 
(ISI)  1 1  /  tm  visibilities  were  best  fit  by  three  spherically  symmet¬ 
ric  concentric  dust  shells.  However,  Monnier  et  al.  (2004)  state 
that  their  Keck-IOTA  data,  which  include  multiple  late-type  stars, 
suggest  that  the  popular  strategy  of  incorporating  multiple  spher¬ 
ically  symmetric  shells  to  model  deviations  from  simple  uniform 
outflow  might  be  misguided,  and  that  dumpiness  and  global  asym¬ 
metry  should  be  considered  more  seriously.  In  fact,  Monnier 
et  al.  (2004)  report  that  the  visibility  data  for  IK  Tau  show  signs 
of  asymmetries  for  the  dust  shell,  but  the  data  are  not  of  sufficient 
quality  to  make  a  definitive  determination. 

Beyond  the  regions  of  the  SiO  masers  and  circumstellar  dust 
shell,  there  is  evidence  for  axisymmetric  structures  in  the  CSE 
of  IK  Tau.  On  large  angular  scales,  Bowers  et  al.  (1989)  found 
that  the  1612  MHz  OH  masers  form  an  elongated  ring  ~3f5  x 
2ff3  with  a  NW-SE  orientation  of  the  long  axis.  The  axisym- 
metry  is  especially  apparent  at  high  velocities  ( V  >  45  km  s_1) 
where  nearly  all  of  the  emission  originates  to  the  SE  and  NW  of 
the  star.  Observations  of  the  22  GHz  H20  masers  with  the  VLA 
(Lane  et  al.  1987;  Bowers  et  al.  1993)  and  the  Multi  Element  Ra¬ 
dio  Linked  Interferometer  Network  (MERLIN;  Yates  &  Cohen 
1994;  Bains  et  al.  2003)  show  similar  axial  symmetry  with  re¬ 
spect  to  the  NW-SE  direction.  VLA  observations  by  Lane  et  al. 
(1987)  show  that  the  low- velocity  masers  have  a  ringlike  distri¬ 
bution,  while  the  high-velocity  masers  have  a  flattened  distribu¬ 
tion  with  a  NW-SE  orientation.  The  overall  size  of  the  shell  was 
roughly  330  mas.  The  VLA  data  of  Bowers  et  al.  (1993)  show 
the  H20  maser  distribution  to  be  ~220  mas  across,  with  a  simi¬ 
lar  elongation  in  the  NW-SE  direction.  The  MERLIN  images  of 
Bains  et  al.  (2003)  show  bright  blue-  and  redshifted  emission 
distributed  over  an  elliptical  region  approximately  1 00  x  200  mas 
superimposed  on  a  fainter  spherical  region  of  emission  from  all 
spectral  channels. 

The  observed  NE-SW  orientation  for  the  elliptical  SiO  maser 
distribution  is  roughly  orthogonal  to  elongations  of  the  H20  and 


OH  maser  distributions.  Bowers  et  al.  (1989)  suggest  that  the 
axisymmetric  structures  seen  in  the  OH  and  H20  maser  shells 
for  certain  late-type  stars  are  a  result  of  the  density  distribu¬ 
tion  of  gas  and  dust  in  the  CSE.  They  conclude  that  the  density 
distribution  is  latitude  dependent,  with  the  total  density  decreas¬ 
ing  as  the  latitude  increases  from  the  rotationally  defined  equa¬ 
torial  plane.  In  stars  in  which  both  the  OH  and  H20  masers  are 
axially  symmetric  (e.g.,  IK  Tau)  Bowers  et  al.  (1989)  suggest 
that  the  distributions  should  be  aligned  parallel  along  the  equa¬ 
torial  plane,  since  the  density  enhancement  is  conducive  to  the 
formation  of  the  H20  masers  (  assuming  these  masers  are  colli- 
sionally  pumped),  and  since  the  OH  masers  are  the  result  of  the 
photodissociation  of  the  H20  by  the  surrounding  UV  radiation. 
Yates  &  Cohen  (1994)  state  that  their  MERLIN  data  are  also 
supportive  of  a  NW-SE  equator  for  IK  Tau  with  the  brightest 
H20  masers  observed  along  this  axis  as  expected.  Assuming 
this  NW-SE  axis  is  parallel  to  the  equator  of  IK  Tau,  as  sug¬ 
gested  by  the  OH  and  H20  maser  observations,  then  the  elon¬ 
gation  we  observe  in  the  SiO  maser  distribution  is  along  an  axis 
roughly  perpendicular  to  the  equator  and  possibly  parallel  to  the 
polar  axis  of  the  star.  The  importance  of  this  geometry  will  be¬ 
come  apparent  in  the  discussion  of  the  velocity  structure  of  the 
SiO  masers  in  the  next  section. 

3.2.  Velocity’  Structure  of  the  SiO  Masers 

The  line-of-sight  (LOS)  velocity  information  derived  from 
our  VLB  A  data  is  shown  in  Figure  3.  The  top  panel  shows  the 
spectrum  of  SiO  maser  emission  from  28  to  42  km  s_1  color 
coded  by  LOS  velocity  in  increments  of  1 .5  km  s  1 .  The  bottom 
panel  of  Figure  3  shows  the  SiO  masers  plotted  in  the  same 
velocity  color  coding  as  the  top  panel,  with  the  color  of  the  ma¬ 
ser  spot  representing  its  corresponding  velocity  in  the  spec¬ 
trum.  It  is  apparent  from  the  figure  that  there  is  a  symmetry 
axis  for  the  velocity  structure  of  the  SiO  maser  emission  toward 
IK  Tau,  with  the  northwestern  side  of  the  shell  dominated  by 
blueshifted  masers  (relative  to  the  assumed  stellar  velocity  of 
35  km  s_1)  and  the  southeastern  side  of  the  shell  dominated  by 
redshifted  masers.  Such  velocity  distributions  in  which  blue- 
and  redshifted  SiO  masers  occur  on  opposite  sides  of  the  star 
have  previously  been  observed  for  NML  Cyg  (Boboltz  &  Marvel 
2000),  R  Aqr  (Hollis  et  al.  2001;  Cotton  et  al.  2004),  and  OH 
231.8+4.2  (Sanchez  Contreras  et  al.  2002).  For  NML  Cyg  and 
R  Aqr,  the  maser  velocity  distributions  suggested  rotation  of  their 
respective  maser  shells.  In  the  case  of  OH  231.8+4.2,  Sanchez 
Contreras  et  al.  (2002)  concluded  that  the  maser  distribution  could 
be  modeled  by  a  rotating,  contracting  torus  of  emission. 

The  ringlike  distributions  typically  observed  in  VLBA  studies 
of  SiO  masers  indicate  that  the  masers  are  amplified  in  a  direction 
perpendicular  to  the  radial  direction  (i.e.,  tangential  amplifica¬ 
tion;  Diamond  et  al.  1994).  This  amplification  occurs  along  the 
longest  path  lengths  parallel  to  the  observer’s  LOS.  Assuming 
this  is  the  case  for  IK  Tau,  then  the  velocity  structure  shown  in 
Figure  3  suggests  a  NW-SE  rotation  of  the  maser  shell  with  the 
axis  of  rotation  near  the  elongation  axis  of  the  distribution.  Such 
a  rotation  axis  is  in  agreement  with  a  NE-SW  equator  for  IK  Tau, 
as  indicated  by  the  previous  OH  and  H20  maser  observations 
mentioned  above. 

Rotation  is  not  a  new  concept  with  regard  to  the  masers  in  the 
CSE  of  IK  Tau.  Bowers  et  al.  (1989)  considered  models  with 
radial  expansion  combined  with  Keplerian  rotation  to  explain 
the  axisymmetric  geometries  observed  for  the  1612  MHz  OH 
masers.  They  concluded  that  although  rotation  could  not  unequiv¬ 
ocally  be  ruled  out,  the  data  did  not  support  rotation  as  a  pri¬ 
mary  component  in  the  kinematics  of  the  OH  shell.  Additionally, 
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Fig.  3. — Structure  of  the  SiO  maser  emission  toward  IK  Tau  as  measured  by 
the  two-dimensional  Gaussian  fit  to  the  masers  in  the  VLBI  channel  maps.  Top 
panel :  Spectrum  formed  by  plotting  component  flux  density  vs.  velocity,  color 
coded  in  1.5  km  s-1  velocity  increments  from  redward  (left)  to  blueward  {right). 
Bottom  panel :  Plots  the  spatial  and  velocity  distribution  of  the  masers.  The  color 
of  each  point  represents  the  corresponding  velocity  bin  in  the  spectrum,  and  the 
size  of  each  point  is  proportional  to  the  logarithm  of  the  flux  density.  Errors  in 
the  positions  of  the  features  are  smaller  than  the  data  points. 

Bowers  et  al.  (1993)  found  that  the  plot  of  radius  versus  veloc¬ 
ity  for  their  H20  maser  data  was  reminiscent  of  models  of  ac¬ 
celerating  or  axisymmetric  outflow  with  a  possible  rotational 
component. 

To  better  characterize  the  velocity  structure  of  the  SiO  emis¬ 
sion,  we  investigated  a  number  of  models  with  velocity  compo¬ 
nents  including  rotation  and  expansion/contraction.  Modeling 
the  three-dimensional  velocity  structure  of  the  SiO  masers  to¬ 
ward  IK  Tau  is  complicated  by  the  elliptical  distribution  of  the 
masers  and  by  a  lack  of  knowledge  of  the  shell  structure  along 
the  observer’s  LOS.  The  elliptical  distribution  of  the  masers 
could  be  modeled  by  either  an  ellipsoidal  shell  or  an  inclined 
torus  similar  to  that  of  Sanchez  Contreras  et  al.  (2002).  We  chose 
an  ellipsoidal  shell,  based  on  the  assumption  that  the  masers 
are  tangentially  amplified  and  the  fact  that  nearly  a  full  ring  of 
emission  is  observed.  In  the  case  of  a  torus  with  tangential  am¬ 
plification,  one  would  expect  to  see  only  two  regions  of  emis¬ 
sion  on  opposite  sides  of  the  star,  as  was  observed  by  Sanchez 
Contreras  et  al.  (2002).  Models  developed  for  IK  Tau  are  based 
on  a  prolate  spheroid  with  the  long  axis  oriented  at  a  position 
angle  of  59°  (as  suggested  by  the  elliptical  fit  shown  in  Fig.  2). 
Since  the  true  three-dimensional  shape  of  the  maser  shell  is 
unknown,  the  axis  of  the  ellipsoidal  model  along  the  observers 
LOS  was  assumed  to  be  the  same  size  as  the  minor  axis  of  the 


model  in  the  plane  of  the  sky.  Additionally,  the  inclination  of  the 
ellipsoid  was  assumed  to  be  zero. 

The  models  were  compared  “by  eye”  to  the  actual  data  to  de¬ 
termine  how  well  they  reproduced  the  true  maser  spatial  and 
velocity  distributions.  To  make  the  comparison  easier,  we  plot¬ 
ted  the  data  and  the  models  in  a  manner  similar  to  that  presented 
in  Hollis  et  al.  (2001),  with  three  different  views.  Figure  4  shows 
the  VLBA  data  plotted  in  the  three  panels  on  the  left.  The  top 
panel  plots  the  maser  distribution  with  the  same  color  coding  as 
that  of  Figure  3,  along  with  a  set  of  axes  (A,  Z)  determined  from 
the  elliptical  fit  of  Figure  2.  The  middle  panel  plots  the  compo¬ 
nent  LOS  velocity  versus  its  position  along  the  equatorial  ( X ) 
axis.  The  bottom  panel  plots  the  component  LOS  velocity  versus 
its  position  along  the  polar  (Z)  axis.  The  right  three  panels  of 
F igure  4  plot  the  equivalent  distributions  derived  from  the  model. 

The  size  and  shape  of  the  model  shell  (i.e.,  major  axis,  minor 
axis,  and  shell  thickness)  were  chosen  to  approximate  the  spatial 
distribution  of  the  SiO  masers  as  seen  in  Figure  3  and  the  top-left 
panel  of  Figure  4.  Maser  velocities  were  first  modeled  with  a 
general  form  of  rotation  given  by 

V(r)ocArp.  (1) 

It  quickly  became  apparent  that  models  including  only  rotation 
do  not  sufficiently  approximate  the  maser  data.  In  the  middle- 
left  panel  of  Figure  4,  we  see  that  there  is  a  significant  variation 
in  position  along  the  A-axis  for  masers  with  small  (<2  km  s~') 
LOS  velocities.  It  is  this  spread  in  position  for  low-velocity 
points  that  is  difficult  to  model  assuming  only  rotation.  For  ex¬ 
ample,  solid-body  rotation  (p  =  1)  yields  a  simple  linear  varia¬ 
tion  of  LOS  velocity  with  position  along  the  A-axis,  regardless 
of  distance  along  the  Z-axis.  Keplerian  rotation  ( p  =  —0.5) 
gives  a  greater  variation  in  A-axis  position  for  higher  velocity 
masers;  however,  at  lower  velocities  the  positions  are  still  con¬ 
fined  to  a  namow  region  in  the  Flos  versus  A  plot  (see  Fig.  3  of 
Hollis  etal.  2001). 

One  way  to  produce  the  observed  spread  in  position  at  small 
LOS  velocities  is  to  add  an  expansion  term  to  the  equation 
describing  the  motion  of  the  masers  in  the  shell.  The  LOS  ve¬ 
locity  for  an  elliptical  shell  that  is  in  rotation  and  also  has  a 
component  of  constant  velocity  expansion,  Vexp,  can  be  written 


V (r)L0S  =  Flpv  +  Y  cos  9cos(j)+  Vsxp  sin  9  cos  <f>,  (2) 

where  the  rotational  term  is  similar  to  that  of  Hollis  et  al.  (2001); 
note  that  the  exponent  q  allows  for  sub-  or  super-Keplerian  ro¬ 
tation.  In  equation  (2),  M  represents  the  mass  of  the  star,  KLpv  is 
the  stellar  velocity,  the  angle  9  is  measured  in  the  equatorial  plane 
of  the  shell,  and  the  angle  (f>  is  the  latitudinal  coordinate  of  the 
shell.  Since  the  mass  of  IK  Tau  is  unknown,  1  M0  was  assumed 
for  the  model.  The  stellar  velocity,  VLFy,  was  assumed  to  be  35  km 
s',  consistent  with  velocities  determined  from  the  midpoints 
of  the  spectra  of  the  1612  MHz  OH  masers,  33.2  km  s  1  ( Bowers 
et  al.  1989),  and  the  H20  masers,  35.4  km  s  1  (Bowers  et  al. 
1993).  Similar  to  the  models  of  Humphreys  et  al.  (1996,  2002), 
maser  gain  length  was  approximated  by  considering  the  distance 
along  the  path  in  which  the  LOS  velocity  changes  by  less  than 
3  times  the  Doppler  width  of  the  line.  Amplification  becomes 
negligible  for  velocity  shifts  greater  than  this  limit  (Humphreys 
et  al.  1996,  2002).  For  SiO  we  estimated  a  Doppler  line  width 
A  Fd  ~  1  km  s  1  for  a  typical  temperature  of  1 500  K  in  the  SiO 
maser  region  ( Doel  et  al.  1 995 ).  The  maser  gain  length  was  further 
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Fig.  4. — Side-by-side  color-coded  plots  of  the  VLBA  data  {left)  and  the  ellipsoidal  shell  model  {right).  Velocity  color  coding  is  the  same  as  for  Fig.  3,  except  that 
the  assumed  stellar  velocity  of  IK  Tau  {Vlpv  —  35  km  s-1)  has  been  subtracted.  The  top  panels  show  the  orientation  of  the  VLBA  data  {left)  and  model  {right) 
distributions  on  the  sky  and  the  assumed  rotational  and  equatorial  axes,  Z  and  X,  respectively,  with  the  rotation  axis  at  P.A.  =  59°.  The  top  left  panel  is  analogous  to 
the  Fig.  3  bottom  panel.  The  middle  and  bottom  panels  show  the  LOS  velocities  plotted  against  the  X  and  Z  axes  of  the  distribution,  respectively,  for  both  the  VLBA 
data  {left)  and  the  model  {right). 


constrained  to  be  greater  than  the  minimum  path  length  of  1 01 3  cm 
necessary  for  SiO  maser  emission  (Elitzur  1982),  and  to  be  less 
than  or  equal  to  the  physical  LOS  path  length  through  the  shell. 

In  practice,  a  random  point  is  chosen  within  the  region  de¬ 
termined  by  the  geometric  parameters  of  the  shell  (i.e.,  size, 
orientation,  and  thickness).  The  LOS  velocity  for  this  point  is 
computed  from  equation  (2).  From  this  point,  distances  forward 
and  backward  along  the  observer’s  LOS  are  computed,  for  which 
the  LOS  velocity  changes  by  1.5  times  the  Doppler  line  width. 
Thus,  we  have  a  line  of  points,  with  the  chosen  random  point  as 


the  center,  for  which  V(r) los  is  within  3  times  the  Doppler  line 
width.  If  both  ends  of  the  line  reside  within  the  maser  shell  and 
the  total  length  of  the  line  is  greater  than  the  minimum  path 
length  required  for  maser  emission,  then  the  center  point  is  kept 
and  is  color  coded  according  to  its  V(r)]  <)S.  Otherwise,  the  point 
is  rejected.  In  this  manner,  the  model  builds  up  a  distribution  of 
points  mapping  both  the  structure  of  the  shell  and  the  maser  ve¬ 
locity  distribution. 

The  model  shown  in  Figure  4  is  meant  to  provide  a  qualita¬ 
tive  match  to  the  IK  Tau  SiO  maser  distribution.  The  parameters 
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TABLE  1 

Ellipsoidal  Shell  Model  Parameters 


Parameter 

Value 

Shell  Geometry: 

Shell  thickness . 

.  2.6  mas 

Position  angle .  59° 

Inclination,  i .  0° 

Distance .  250  pc 

Stellar  mass,  M .  1  M0 

Stellar  velocity,  Flpv .  35.0  km  s-1 

Minimum  maser  path  length .  lx  1013  cm 

Doppler  line  width,  AFd .  1.0  km  s-1 

Expansion  velocity,  Vexp .  3.5  km  s-1 

Exponent,  <7 .  1.09 


for  the  model  are  presented  in  Table  1 .  Inner  distances  of  1 6 
and  28  mas  were  used  for  the  semiminor,  A0,  and  semimajor,  B0, 
axes,  respectively,  with  a  shell  thickness  of  3  mas.  As  men¬ 
tioned  previously,  the  position  angle  for  the  ellipse  was  set  to 
59°.  The  model  is  not  particularly  sensitive  to  this  angle,  and 
changes  of  515°  do  not  alter  the  results  significantly.  The  two 
parameters  that  affect  the  velocity  structure  are  the  exponent 
q  and  the  expansion  velocity  Vexp .  Through  an  iterative  process 
of  adjusting  parameters  and  qualitatively  comparing  the  model 
to  the  data,  values  of  q  =  1.09  and  Vexp  =  3.5  km  s_1  were  de¬ 
termined.  Coincidentally,  the  value  for  q  is  equivalent  to  that 
found  for  R  Aqr  by  Hollis  et  al.  (2001).  As  in  the  model  for 
R  Aqr,  the  model  for  IK  Tau  is  sensitive  to  small  (^0.01)  var¬ 
iations  in  the  value  of  q.  Unlike  the  R  Aqr  case,  the  mass  of 
IK  Tau  is  not  well  constrained,  and  variations  in  the  assumed 
value  forMcould  allow  for  slight  variations  in  the  value  of  q.  To 
produce  a  model  similar  to  that  shown  in  Figure  4  with  pure 
Keplerian  motion  (q  =  1 )  would  require  a  stellar  mass  of  M  k 
0.1  Mq.  Similarly,  the  model  is  sensitive  to  the  expansion  veloc¬ 
ity  Eexp.  Changes  on  the  order  of  ~1.0  km  s_1  produce  signif¬ 
icant  changes  in  the  velocity  distributions  of  the  masers. 

Additionally,  there  is  a  trade-off  between  the  model  expan¬ 
sion  velocity  and  the  quenching  of  the  masers  caused  by  the 
Doppler  line  width  limitation.  Larger  values  for  Kexp  tend  to 
spread  the  distribution  in  the  middle  panel  of  the  model  plot,  but 
only  if  A  VD  is  increased  to  allow  the  greater  change  in  velocity 
and  resulting  longer  gain  lengths.  The  longer  path  lengths,  in 
turn,  create  more  masers  in  front  of  and  behind  the  star,  a  char¬ 
acteristic  that  is  clearly  not  evident  in  these  or  any  other  SiO  ma¬ 
ser  observations  of  which  we  are  aware.  As  the  top  right  panel  of 
Figure  4  shows,  the  density  distribution  of  masers  resulting  from 
the  path  length  considerations  is  characterized  by  tangentially  am¬ 
plified  masers  along  the  limb  of  the  shell.  In  order  to  limit  maser 
emission  along  the  LOS  to  the  star,  the  maser  shell  thickness 
must  be  less  than  or  approximately  equal  to  the  maser  minimum 
path  length.  Since  we  assumed  a  minimum  maser  path  length  of 
1013  cm  (~2.7  mas  at  the  assumed  distance  to  IK  Tau),  the  shell 
thickness  was  set  to  2.6  mas.  Without  increasing  the  minimum 
required  path  length,  increasing  the  maser  shell  thickness  will  pro¬ 
duce  masers  along  the  LOS  to  the  star,  thus  filling  in  the  distribu¬ 
tion  in  the  top-right  panel  of  Figure  4. 

The  expansion  velocity  Vexp  =  3.5  km  s_1  used  in  the  model 
is  consistent  with  observational  evidence  as  found  in  proper- 
motion  studies  of  SiO  masers.  Boboltz  et  al.  (1997)  measured 
the  proper  motions  of  the  SiO  masers  toward  R  Aqr  and  found 


an  overall  infall  of  the  masers  (shell  contraction)  with  a  velocity 
of  ^4  km  s_I  over  stellar  phases  <j>  =  0.78-1.04.  Diamond  & 
Kemball  (2003),  in  an  extensive  study  of  the  Mira  variable  TX 
Cam,  determined  the  dominant  mode  to  be  expansion  of  the  SiO 
masers,  with  radial  proper-motion  magnitudes  of  ~5  -1 0  km  s~ 1 
between  stellar  phases  </>  s=s  0.5 -1.5. 

From  our  single  epoch  of  SiO  maser  observations,  it  is  dif¬ 
ficult  to  tell  whether  the  model  expansion  term  represents  a 
physical  expansion  or  contraction  of  the  maser  shell.  This  is  be¬ 
cause  both  expansion  and  contraction  produce  nearly  identi¬ 
cal  distributions  in  the  both  the  middle-  and  lower-right  panels 
of  Figure  4.  The  only  difference  between  a  positive  or  negative 
expansion  term  is  observable  in  the  top-right  panel  of  Figure  4. 
For  example,  if  we  consider  only  rotation  of  the  shell,  we  find 
that  all  of  the  redshifted  masers  occur  along  the  positive  X-axis 
(the  SE  side  of  the  shell ),  and  all  of  the  blueshifted  masers  occur 
along  the  negative  X-axis  (NW  side  of  the  shell),  and  the  ve¬ 
locities  are  symmetric  about  the  plane  of  the  sky.  The  addition 
of  a  positive  expansion  term  (as  was  done  for  our  model )  serves 
to  rotate  the  blueshifted  velocities  out  of  the  plane  of  the  sky 
toward  the  front  of  the  shell  closest  to  the  observer  and  the  red- 
shifted  velocities  away  from  the  observer.  However,  the  actual 
masers  are  still  confined  to  a  region  around  the  plane  of  the  sky; 
thus,  the  shift  in  the  velocity  structure  is  most  easily  observed 
near  the  poles  along  the  Z-axis  where  we  see  in  the  top-right 
panel  of  Figure  4  that  some  low-velocity  blueshifted  masers 
(green)  have  been  shifted  toward  the  positive  X-axis.  Contrac¬ 
tion  of  the  shell  has  the  opposite  effect,  rotating  the  redshifted 
velocity  structure  out  of  the  plane  toward  the  front  of  the  shell 
and  the  blueshifted  velocities  toward  the  back  of  the  shell.  It  is 
practically  impossible  to  tell  from  our  maser  data  whether  ex¬ 
pansion  or  contraction  is  the  case  for  the  IK  Tau  shell.  There 
seems  to  be  a  slight  preference  for  blueshifted  masers  in  general 
and  a  tendency  for  some  blueshifted  masers  along  the  positive 
X-axis  in  the  top-left  panel  of  Figure  4,  but  this  is  highly  de¬ 
pendent  on  the  assumed  stellar  velocity  of  35  km  s_1  and  the 
choice  of  position  angle  for  the  major  axis  of  the  ellipse. 

The  maser  shell  rotational  velocities  produced  in  the  model 
are  not  considered  extraordinary.  The  equatorial  velocity  at 
the  distance  of  the  mean  semiminor  axis  of  16  mas  is  ~3.6  km 
s_1  with  a  corresponding  rotational  period  of  ~33  yr.  At  the 
UD  photospheric  distance  of  10. 1  mas,  the  equatorial  velocity  is 
^4.6  km  s_1  with  a  corresponding  rotational  period  of  ~16  yr. 
These  rotational  velocities  are  consistent  with  observational  mea¬ 
surements  for  R  Aqr  (~3-5  km  s  1 ;  Hollis  et  al.  2001;  Cotton 
et  al.  2004)  and  slightly  lower  than  measurements  for  NML  Cyg 
(~  1 1  km  s“ 1 ;  Boboltz  &  Marvel  2000)  and  OH  23 1 . 8+4.2  (-7- 
10  km  s_1;  Sanchez  Contreras  et  al.  2002).  Significant  stellar 
rotation  velocities  (10-100  s  of  km  s_1)  during  the  asymptotic 
giant  branch  (AGB)  phase  have  been  invoked  in  theoretical  stud¬ 
ies  of  chemical  mixing  and  the  nucleosynthesis  of  heavy  elements 
in  convective  envelopes  (Langer  et  al.  1999;  Herwig  et  al.  2003; 
Siess  et  al.  2004).  Stellar  rotation  is  also  a  factor  in  models  de¬ 
scribing  the  shaping  of  bipolar  and  elliptical  planetary  nebulae, 
and  Garcia-Segura  et  al.  (1999)  propose  that  AGB  stars  above 
~1.3  Mq  can  spin  up  their  extended  circumstellar  envelopes  to 
rotational  speeds  on  the  order  of  a  km  s  1 .  Stellar  rotation  ve¬ 
locities  of  ~3— 7  km  s  1  are  used  in  such  models,  which  result  in 
equatorial  density  enhancements  and  eventually  bipolar  nebulae 
(Garcia-Segura  et  al.  1999).  Interestingly,  the  collective  maser 
observations  of  IK  Tau,  which  suggest  rotation,  axial  symmetry, 
and  equatorial  density  enhancement,  may  be  indicating  that  IK 
Tau  is  on  its  way  to  becoming  such  a  bipolar  planetary  nebula. 
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4.  SUMMARY 

We  have  imaged  the  v  =  1,  J  =  1-0  SiO  maser  emission 
toward  the  Mira  variable  IK  Tau.  We  find  that  the  SiO  masers 
are  distributed  in  an  elongated  ring  roughly  58x32  mas  (15x 
8  AU  at  the  adopted  distance  of  250  pc)  with  the  major  axis  of 
the  ellipse  at  a  position  angle  of  ^59°.  The  NE-SW  elongation 
of  the  SiO  maser  distribution  is  roughly  perpendicular  to  the 
NW-SE  symmetry  axis  previously  determined  for  the  H20  and 
OH  maser  distributions.  This  elongation  axis  also  provides  an 
axis  of  symmetry  for  the  LOS  velocity  structure  of  the  SiO  ma¬ 
sers.  Masers  having  velocities  blueshifted  with  respect  to  the 
assumed  center  velocity  of  35  km  s_1  were  found  to  lie  on  the 
northwestern  side  of  the  shell,  while  the  southeastern  side  of 
the  shell  is  dominated  by  redshifted  SiO  masers.  Under  the  as¬ 
sumption  of  tangentially  amplified  masers,  this  velocity  dis¬ 
tribution  implies  a  NW  to  SE  rotation  of  the  SiO  shell  about 
a  NE-SW  rotation/polar  axis.  Such  a  polar  axis  would  be  in 
agreement  with  models  invoked  to  explain  the  H20  and  OH 
maser  geometries,  which  suggest  a  NW-SE  equator  for  IK  Tau. 
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The  inherent  transient  and  irregular  nature  of  SiO  maser 
emission  in  the  atmospheres  of  long-period  variables  limits  the 
comparison  of  models  to  a  single  epoch  of  VLBI  data.  How¬ 
ever,  we  have  attempted  to  approximate  the  IK  Tau  observa¬ 
tions  presented  here  assuming  an  ellipsoidal  model  of  the  SiO 
maser  shell.  We  favor  a  model  that  includes  both  rotation  of  the 
form  V  cx  1/r109  and  constant-velocity  expansion/contraction 
as  a  representation  of  the  velocity  field  in  the  CSE  of  IK  Tau 
at  the  distance  of  the  SiO  maser  shell.  This  model  appears  to 
reproduce  the  global  characteristics  of  the  observed  SiO  masers 
in  terms  of  structure  and  velocity  distribution.  A  multiepoch 
VLBI  monitoring  program  to  determine  the  SiO  maser  proper 
motions  could  verify  this  scenario  and  would  greatly  clarify  our 
understanding  of  the  CSE  of  IK  Tau. 
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